Timescales of Turbulent Relative Dispersion 
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Tracers in a turbulent flow separate according to the celebrated t 3 ^ 2 Richardson-Obukhov law, 
which is usually explained by a scale-dependent effective diffusivity. Here, supported by state-of- 
the-art numerics, we revisit this argument. The Lagrangian correlation time of velocity differences 
is found to increase too quickly for validating this approach, but acceleration differences decorrelate 
on dissipative timescales. This results in an asymptotic diffusion oc i 1//2 of velocity differences, so 
that the long-time behavior of distances is that of the integral of Brownian motion. The time of 
convergence to this regime is shown to be that of deviations from Batchelor's initial ballistic regime, 
given by a scale-dependent energy dissipation time rather than the usual turnover time. It is finally 
argued that the fluid flow intermittency should not affect this long-time behavior of relative motion. 



Turbulence has the feature of strongly enhancing the dis- 
persion and mixing of the species it transports. It is 
known since the work of Richardson [l| that tracer parti- 
cles separate in an explosive manner oc t 3 ^ 2 that is much 
faster and less predictable than in any chaotic system. 
While little doubt remains about its validity in three- 
dimensional homogeneous isotropic turbulence, observa- 
tions of this law in numerics and experiments are dif- 
ficult, as they require a huge scale separation between 
the dissipative lengths, the initial separation of tracers, 
the observation range and the integral scale of the flow 
0, Q ■ Much effort has been devoted to test the univer- 
sality of this law, which was actually retrieved in various 
turbulent settings, such as the two-dimensional inverse 
cascade buoyancy- driven flows 5], and magneto- 
hydrodynamics [6] . At the same time, breakthroughs on 
transport by time-uncorrelated scale- invariant flows have 
strenghtened the original idea of Richardson that this 
law originates from the diffusion of tracer separation in a 
scale-dependent environment 0- As a result, the physi- 
cal mechanisms leading to Richardson-Obukhov t 3 / 2 law 
are still rather poorly understood and many questions re- 
main open on the nature of subleading terms, the rate of 
convergence and on the effects of the intermittent nature 
of turbulent velocity fluctuations d, Q . 

Turbulent relative dispersion consists in understand- 
ing the evolution of the separation 6x(t)=Xi(t) — X^it) 
between two tracers. Richardson's argument can be 
reinterpreted by assuming that the velocity difference 
5u(t) = u(Xi, t) — u(X%, t) has a short correlation time. 
This means that the central-limit theorem applies and 
that, for sufficiently large timescales, 



dSx 



(1) 



where £ is the standard three-dimensional white noise, 
U T U = {Su <g> Su) the Eulerian velocity difference corre- 
lation tensor, and tl the Lagrangian correlation time of 
velocity differences between pair separated by &c=|&c|. 
As stressed by Obukhov , when assuming Kolmogorov 
1941 scaling, tl ~ 8x 2 / 3 , U ~ Sx 1 / 3 , and the Fokker- 



Planck equation associated to (p} exactly corresponds to 
that derived by Richardson for the probability density 
p(5x,t). It predicts in particular that the squared dis- 
tance (\5x{t)\ 2 ) ro averaged over all pairs that are initially 
at a distance |&c(0)| — ro has a long-time behavior oc t 3 
that is independent on ro . This memory lost on the initial 
separation can only occur on time scales longer than the 
correlation time tt, (ro) ~ r 2 / 3 of the initial velocity differ- 
ence. For times t -C Tt(ro), one cannot make use of the 
approximation (p} as the velocity difference almost keeps 
its initial value. This corresponds to the ballistic regime 
(|fe(<)-^(0)| 2 ) ro ~t 2 S 2 {r ), where S 2 (r) = (|<5«| 2 ) ^ 
the Eulerian second-order structure function over a sep- 
aration r, introduced by Batchelor [ll|. The diffusive 
approach JT]) can however be modified to account for the 
ballisitic regime [12j |. Nevertheless a short-time correla- 
tion of velocity differences can hardly been derived from 
first principles and seems to contradict turbulence phe- 
nomenology. Indeed, as stressed in [7], if Sx grows like 
t 3 / 2 , the Lagrangian correlation time tt, is of the order of 
8x 2 / 3 ~t, so that the velocity difference correlation time 
is always of the order of the observation time. Despite 
such apparent contradictions, Richardson diffusive ap- 
proach might be relevant to describe some intermediate 
regime valid for large-enough times and typical separa- 
tions. Several measurements show that the separations 
distribute with a probability that is fairly close to that 
obtained from an eddy-diffusivity approach 0, [l3, 14 1. 

To clarify when and where Richardson's approach 
might be valid, it is important to understand the 
timescale of convergence to the explosive t 3 law. Much 
work has recently been devoted to this issue: it was for in- 
stance proposed to make use of fractional diffusion with 
memory [1 51 ] . to introduce random delay times of con- 
vergence to Richardson scaling [16], or to estimate the 
influence of extreme events in particle separation [l7j ]. 
All these approaches consider as granted that the final 
behavior of separations is diffusive. As we will see here, 
many aspects of the convergence to Richardson's law for 
pair dispersion can be clarified in terms of a diffusive 
behavior of velocity differences. 
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To address such issues, we make use of direct numerical 
simulations. For this, the Navier-Stokes equation with a 
large-scale-forcing is integrated in a periodic domain us- 
ing a massively parallel spectral solver at two different 
resolutions. Table fl] summarizes the parameters of the 
simulations (see [18| for more details). In each case, the 
flow is seeded with 10 7 Lagrangian tracers. Their posi- 
tions, velocities, and accelerations are then stored with 
enough frequency to study relative motion. 
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TABLE I: Parameters of the numerical simulations. N is the 
number of grid points, Rx the Taylor-based Reynolds num- 
ber, v the kinematic viscosity, e the averaged energy dissipa- 
tion rate, u rms the root-mean square velocity, r) = (u'^/e) 1 ^ 4 
the Kolmogorov dissipative scale, t v = (u/e) 1 ^ 2 the associated 
turnover time, L = u 3 ma /e the integral scale and T = L/u vulB 
the associated large-scale turnover time. 
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FIG. 1: Time-evolution of the mean-square separation for 
Rx — 730 and different initial separations. The dashed line 
represents the two leading terms of the ballisitic behavior 
@. The solid line is a fit to the Richardson regime Q 
with g — 0.52 and C — 1.6. Inset: to as a function of ro in 
dissipative-scale units. The solid line is an Eulerian average, 
the circles are Lagrangian measurements and the dashed line 
is the turnover time r(ro). 



We first report results on the behavior of the separa- 
tion Sx(t) as a function of time. Following [l3j], a Taylor 
expansion at short times leads to 

(\Sx(t)-6x(0)\ 2 ) ro = t 2 5 2 (r )+t 3 (Su • (5a)+0(r 4 ) , (2) 

where S2 (r) = (|<5it| 2 ) is the second-order structure func- 
tion, (•) denote Eulerian averages, and Sa(t) = a(Xi,t) — 
a(X2,t) is the difference of the fluid acceleration sampled 
by the two tracers (where the notation a = dtu+u ■ Vit). 
As long as the term oc i 2 is dominant, the tracers sepa- 
rate ballistically. Expansion © clearly fails for t ps to = 
5a(n))/| (Su -Sa)\. It is known @, [l9| that for separa- 
tions in the inertial range (Su ■ 8a) = — 2e, which is noth- 
ing but a Lagrangian version of the 4/5 law. This implies 
that the ballistic regime ends up at times of the order of 



to = S 2 (r )/(2e). 



(3) 



This timescale can be interpreted as the time required 
to dissipate the kinetic energy contained at the scale ro- 
We thus expect it to be equal to the correlation time of 
the initial velocity difference, to differs from the turnover 
time r(r ) = ^/[^(ro)] 1 / 2 defined as the ratio between 
the separation ro and the typical turbulent velocity at 
that scale. When Kolmogorov 1941 scaling is assumed, 
these two time scales have the same dependency on ro. 
However, usual estimates of the Kolmogorov constant 
lead to to/r(ro) ~ 20. Also, note that intermittency cor- 
rections to the scaling behavior of S2 should in principle 
decrease this ratio. Figure[T]represents the mean-squared 
displacement rescaled by ^5a(ro) as a function of t/to, 
for various values of the initial separation ro. In such 
units and when ro is far in the inertial range, all mea- 
surements collapse onto a single curve. The subleading 
term oc t 3 in @ is relevant for times t < 0.01 to- 



The data collapse extends to times larger than to when 
the mean squared separation tends to Richardson t 3 
regime. This unexpected fact implies that to is not only 
the timescale of departure from the ballistic regime, but 
also that of convergence to Richardson's law. More pre- 
cisely, numerical data suggest that for t 3> to 

(|^(t)-^(0)| 2 > rD =.get 3 [l + C(t /t)]+h.o.t.. (4) 

The constant C does not strongly depend on the 
Reynolds number. Systematic measurements as a func- 
tion of the initial separation show that C is negative 
when ro is of the order of the Kolmogorov scale r/. The 
convergence to Richardson law is then from below and 
is thus contaminated by tracer pairs which spend long 
times close together before sampling the inertial range; 
this is consistent with the findings of 17(. When ro is 
far-enough in the inertial range, C ~ 1.6 becomes inde- 
pendent on the initial separation and the convergence to 
Richardson law is from above. One finds that C = for 
ro ~ 4r7; the only subleading terms in (0| are then of 
higher order, so that the mean-squared separation con- 
verges faster to Richardson regime. Such an initial sep- 
aration could be an "optimal choice" to observe the t 3 
behavior in experimental settings. 

To understand why the timescale of convergence to 
Richardson law is of the order of to, let us examine the 
timescales entering the relative dispersion process. As al- 
ready stated, the velocity difference Su between the two 
tracers stays correlated over a time that increases too fast 
with the separation, making difficult to justify the diffu- 
sive approach (jTJ) . However, it is known that turbulent 
acceleration, which is a small-scale quantity, is correlated 
over times that are of the order of the Kolmogorov 
turnover time [20j |. Its amplitude is rather correlated on 
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times of the order of the forcing correlation time, but 
this does not alter the argument below. Figure [2] repre- 




FIG. 2: Lagrangian time autocorrelation of the acceleration 
difference 5a for various ro and R\ — 730. Inset: for the 
same separations ro, variance of the acceleration difference 
amplitude conditioned on the longitudinal velocity difference 
5«" as a function of the local dissipation rate [<ki"] 3 /ro- 



sents the Lagrangian autocorrelation of the difference of 
acceleration 5 a between two tracers. We clearly see that 
the components of this quantity decorrelate on times of 
the order of t v . This suggests applying the central- limit 
theorem, so that for separations in the inertial range and 
on timescales much longer than the t v , the difference of 
acceleration between two tracers can be approximated by 
a delta-correlated-in-time random process. We thus have 



dSx 
~dT 



= 8u, with — — = 8a~ x /T}° c A(5x, 5u) £(t), (5) 



where A is defined as A T A = (5a (g) 5a\5x,5u), £ is 
the three-dimensional white noise, and the product is 
here understood in the Stratonovich sense. The idea 
of assuming uncorrelated accelerations is common to 
many stochastic models for turbulent dispersion (see, 
e -g- ; [1, 21[). However, this model does not require any 
Eulerian input and involves a multiplicative noise (A de- 
pends on 5u). Dimensional arguments indicate that the 
local Kolmogorov time r^ oc and the acceleration ampli- 
tude A = \A\ depends only on the viscosity v and on 
the local energy dissipation rate ei oc . We thus have 



T \o C _ ^/ 2 e lo l /2 and A ~ ^ 1/4 ef^. These estimates 
predict that the multiplicative term in ([5]) behaves as 
[t^ 00 ] 1 / 2 ^ ~ e^ 2 . Interestingly this quantity is indepen- 
dent on v and is thus expected to have a finite limit 
at infinite Reynolds numbers. Phenomenological argu- 
ments suggest that for typical values of the velocity dif- 
ference 5u, the local dissipation rate can be written as 
eioc ~ [<5u"] 3 /<5x, where 5u" = Sx ■ 5u/5x is the lon- 
gitudinal velocity difference between the tracers. When 
<5?i" = 0, the local dissipation rate does not vanish but can 
be estimated through an averaged contribution of larger 



eddies, leading to ei oc ~ e, the averaged energy dissipa- 
tion rate. These estimations have been tested against 
numerical simulations: the inset of Fig. [2] shows the vari- 
ance of the acceleration differences conditioned on the 
longitudinal velocity difference for various separations. 
Up to some statistical errors, it seems that data are in 
rather good agreement with the phenomenological pre- 
diction which is shown as a dashed line. Finally such 
dimensional considerations lead to model the large-time 
evolution of tracer separation as 



&5x 
~~dT 



= 5uK 



dSi 



dt 



e + a ■ 



[fa"] 3 
8x 



1/2 



£(*), (6) 



where a is a positive parameter. Again here the multi- 
plicative noise is understood with Stratonovich conven- 
tion. When rewriting it in the Ito sense, the additional 
drift that appears introduces a "correlation time" equal 
to the instantaneous turnover time Sx/Svfi. Preliminary 
studies of (j6]) showed that its solutions follow a ballistic 
regime at short times and behave according to Richard- 
son law, i.e. (5x 2 ) ~ i 3 at large times. In this stochastic 
model, the local dissipation [<5uII] 3 /<5:e tends to a con- 
stant at large times, so that in the asymptotic regime, 
the velocity difference obeys an equation of the form 
d<5it"/di oc £(t) and thus diffuses. So far, we have only in- 
vestigated the one-dimensional version ([6]) of the model. 
Extension to higher dimensions requires accounting for 
incompressibility and is the subject of ongoing work. 
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FIG. 3: Time evolution of the averaged longitudinal velocity 
for R\ — 730 and different ro (same symbols as in Fig. [TJ. 
The short-time prediction (0 is shown as a dashed line. The 
diffusive behavior (\5u\ 2 ) ro ~ S^ro) + 2.3 e t is represented as 
a dash-dotted line. Inset: time evolution of ([5u^] 3 /5x) ra ; the 
dashed line corresponds to the value 6e. 



To address the relevance of such a model to real flows, 
we turn back to the analysis of simulation data. Figure [3] 
shows the time evolution of (\5u(t)\ 2 ) ro for various values 
of ro. At small times this quantity slightly decreases 
because the subleading term is negative. We indeed have 
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5u(t) ~ 8u(0) +t8a(0), so that the ballistic regime reads 

(\6u{t)\ 2 ) ra =5 a (r )(l-2t/to) + h.o.t. (7) 

Again, the subleading terms are relevat for times t < 
0.01 to- Figure [3] also shows that at large times the mean- 
squared velocity difference looses dependence on ro and 
grows oc et. In addition, as seen from the inset of Fig. [31 
the averaged local dissipation rate ([Su^] 3 /5x) ro along 
particle pairs approaches a positive constant ~ 6 e (inde- 
pendently on R\) on times of the order of t v . This con- 
firms the relevance of the mechanisms described above in 
terms of a stochastic equation for the velocity differences. 

Numerical results indicate that the time to controls the 
convergence to a diffusive regime for initial separations ro 
far enough in the inertial range. This can be explained 
by the following argument. As ([Su^] 3 /Sx) ro becomes 
constant on a short timescale, one expects that 

(\Su(t)\ 2 ) ro ~S 2 (r ) + Det for t » t„, (8) 

where D is a positive constant (for both Reynolds num- 
bers, we observe D s» 2.1). By balancing the diffusive 
term with the initial mean-squared velocity difference 
(\Su(0)\ 2 ) ra = S^ro), we find again that the former is 
dominant for times t much larger than to- The diffu- 
sive behavior of velocity differences is thus reached at 
times of the order of to and this explains in turn why this 
timescale is that of convergence to Richardson's regime. 

Let us summarize here our findings. In this work we 
give some evidence that the Richardson explosive regime 
(|(5a;| 2 ) oc t 3 for the separation between two tracers in 
a turbulent flow originates from a diffusive behavior of 
their velocity difference rather than from dimensional ar- 
guments or equivalently a scale-dependent eddy diffusiv- 
ity for their distance. This leads on to reinterpret the 
t 3 law as that of the integral of Brownian motion. Such 
an argument is supported by two observations. First, the 
acceleration difference has a short correlation time (of the 
order of the Kolmogorov dissipative timescale) and can 
be approximated as a white noise. Second, the ampli- 
tude of this noise solely depends on the local dissipation 
rate ([5u^] 3 /5x) ro , which becomes constant also on short 
timescales. These considerations allow us to show that 
the time to of convergence to Richardson's law is equal to 
that of deviations from Batchelor's ballistic regime. This 
time, which reads to = 52(r )/(2e), is the time required 
to dissipate the kinetic energy contained at a scale equal 
to the initial separation between tracers. 

The interpretation of Richardson's law as the diffusion 
of velocity differences strongly questions possible effects 
of fluid-flow intcrmittcncy on trajectory separation. In- 
deed, considerations on velocity scaling, which are pri- 
mordial in approaches based on eddy diffusivity, are ab- 
sent from the arguments leading to a diffusive behavior 
of 5u. Hence, we expect the separation Sx to follow a 
self-similar evolution in time, independently on the or- 
der of the statistics. Intermittency will however affect s 



directly the time of convergence to such a regime. More 
frequent violent events (of tracer pairs approaching or 
fleeing away in an anomalously strong manner) will result 
in longer times for being absorbed by the average. Such 
arguments do not rule out the possibility of having inter- 
mittency corrections when interested in other observables 
than moments of the separation, as it is for instance the 
case for exit times [8(. Such issues will certainly gain 
much from a systematic study of multi-dimensional gen- 
eralizations of the stochastic model introduced here. 
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